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We show that the resistive anisotropy of an anisotropic medium is determined by the ratio of the phase coherence 
lengths. In layered crystals in which the interlayer transport is incoherent, the out-of-plane phase coherence length 
is fixed and temperature independent. This leads to a temperature dependent resistive anisotropy and to the 
coexistence of metallic in-plane and non-metallic out-of-plane conductivities. Our approach provides a description 
of the c-axis conductivity in the highly nonclassical regime, characteristic of layered cuprates. 
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The normal state of highly anisotropic crystals such 
as high-Tc cuprates exhibits a number of unusual fea- 
tures. The coexistence of metallic in-plane (p a ) and non- 
metallic out-of-plane (p c ) resistivities within a range of 
doping is one. A crossover from metallic to non- 

metallic temperature (T) dependence is another [|],[|. 
From a theoretical point of view, there is a great deal of 
interest in exploring the consequences of the idea that in- 
terlayer transport in these crystals is incoherent, so that 
successive interlayer transitions are uncorrelated as evi- 
denced by the analysis of dc and infrared conductivities 



The coherence of the wave function of the charge car- 
riers can be described in terms of the phase coherence 
length, hereafter called the Thouless length (TL), the 
distance electrons travel between dephasing inelastic col- 
lisions ||. An anisotropic medium, naturally, is char- 
acterized by an anisotropic TL. Below we show that, in 
general, the anisotropy a a /a c is determined by the ra- 
tio of TL in the respective directions. The incoherence of 
the out-of-plane transport in layered crystals implies that 
TL in the c-direction is temperature independent, equal 
to the spacing between neighboring layers. As a result, 
the temperature dependence of p c j p a is completely de- 
termined by the T-dependence of the in-plane TL. This 
can lead to coexistence of metallic p a and non-metallic 
p c - We also explore the consequences of the idea that the 
evolution of the conductivities (a a and er c ) with temper- 
ature and doping can be described through their depen- 
dence on the in-plane Thouless length. One of the most 
interesting conclusions that can be drawn from this anal- 
ysis is that there may exist a unifying description of the 
conductivity of layered crystals at different doping levels. 

Let us consider an anisotropic medium characterized 
by three different TL: £ x ,£ y , and £ z . The phase-coherent 
volume of such a medium, which may be considered as a 
block with sizes equal to the respective TL (hyperblock) 
has isotropic conductance; i.e., g x — g y — g z = 9- One 
way to prove this is by using the original Thouless idea 
that the conductance of a block is 



where i = {x, y, z), dN/dE is the total number of states 
inside the block per unit energy, and (AE)i is the mean 
fluctuation in energy levels caused by replacing periodic 
by antiperiodic boundary conditions in the direction of 
the current. The value of (AE)i is determined by the 
sensitivity of the energy levels to the boundary condi- 
tions, which in turn depends on the phase memory. Since 
the decoherence of the wavefunction occurs within the 
boundaries of the hyperblock, (AE)i should be isotropic. 
Also, (AE)i is determined by the time it takes for a par- 
ticle to cross the system in a given direction ||. By 
definition, it takes an electron the same dephasing time 
t v to cross the hyperblock in each direction, so that 
(AE)i ~ H/t v . 

A macroscopic block {L x , L y , L z } obtained by fitting 
together N 3 hyperblocks (L x /£ x — L y /£ y — L z /t z = 
N 3> 1) also has isotropic conductance G s» Ng which 
can be expressed in terms of the components of the con- 
ductivity tensor: G = a x L y L z / L x = a y L x L z / L y — 
a z L y L x / L z . As a result, we arrive to the following rela- 
tionship between conductivities: 
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This is a general result, valid for all media, not only for 
layered crystals. For example, in an anisotropic Fermi liq- 
uid li oc Vf^Tip, where Vpj is a component of the Fermi 
velocity. Then, Eq. (2) reduces to the ratio of the ef- 
fective masses. In the case of anisotropic diffusion, the 
Thouless lengths £f oc DiT v , where Di is a component 
of the diffusion tensor. For such a system, Eq. (2) is 
equivalent to another known result <Xj / aj = Di /Dj . 
Finally, for the case of variable range hopping (VRH), 
Eq. (2) was obtained through a different and less general 
approach, where is the average hopping distance in the 
respective direction [^0| . 

We apply Eq. (2) to highly anisotropic layered crystals 
in which electrons lose coherence in the c-direction over 
the smallest possible distance, the interlayer spacing £q 
(unidirectional Ioffe-Regel limit, £ z — const — £q). For 
simplicity, the planes are considered to be isotropic; i.e., 
£ x = l y = £. Under these conditions, Eq. (2) gives 
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Thus, the out-of-plane conductivity, which has been one 
of the most enigmatic features of the normal state of 
cuprates, is completely determined by the in-plane con- 
ductivity and the in-plane Thouless length. 

Typically, quantum phenomena (e.g., superfluidity 
and superconductivity) reveal themselves macroscopi- 
cally through phase coherence established over macro- 
scopic distances. Here, anomalously strong interlayer 
decoherence in the normal state (of still unknown ori- 
gin) has macroscopic consequences such as a strongly 
T-dependent resistive anisotropy which reflects the T- 
dependence of the in-plane phase coherence length. 
Equation (3) cannot be obtained from the quasiclassi- 
cal kinetic equation because c-axis transport is highly 
nonclassical even when the in-plane transport can be, in 
principle, treated quasiclassically. 

As a first example of the applicability of Eq. (3) to 
layered crystals, we show below that for optimally and 
nearly optimally doped cuprates, Eq. (3) along with dif- 
fusive in-plane transport gives a description of the out- 
of-plane resistivity which agrees very well with experi- 
ment and explains the observed correlation between p a 
and p c . While Eqs. (2) and (3) are fundamental, the 
T-dependence of the in-plane TL is not, since TL can be 
determined by a variety of dephasing processes with one 
of them dominant in a certain T range for a given system. 
Here we assume that the strongest dephasing process is 
due to the thermal spread of energies, so that TL is deter- 
mined by the thermal diffusion length. [ Support for this 
assumption is provided by Ref. [fill which shows that the 
electron-electron interaction may be equivalent to the in- 
teraction of the electrons with the thermal fluctuations of 
the electromagnetic waves.] Then, the in-plane TL may 
be written as 
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Here we have added the empirical cutoff length £ to 
the conventional definition because in a crystal, even at 
high temperatures, the phase coherence length cannot be 
shorter than a certain finite length, e.g. interatomic dis- 
tance. Since elastic collisions do not lead to the loss of 
coherence, more realistically £ is comparable to the elas- 
tic mean free path [at least far enough from the metal- 
insulator transition (MIT) as we discuss later]. 

With the phenomenological expression for the in-plane 
resistivity, p a = (3 a + a a T H, Eqs. (3) and (4) give: 
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and To = hD/i^ks- Two important aspects are worth 
noting. First, a large constant p c appears even in the 
absence of the intercept (3 a (if p a — a a T, then p c — 
P c + a c T). This apparent "residual resistivity" violates 
Matthiasen's rule, because it does not scale to zero with 
the concentration of impurities. Second, a finite positive 
intercept f3 a in p a translates into a non-metallic term 
7c/T, so that p c has a minimum. Such correlations be- 
tween p a and p c exist in optimally and slightly under- 
doped cuprates fll. From the data of Y Ba-iCuzOi-s 
p"2[ for which a a Rj 0.5 pfl cm/K, P c rj 1.1 mfl cm, 
and a c RJ 12.5 pfl cm/K, we estimate To RJ 2200 K, 
D rj 4 cm 2 Is, and t\ = i Q {a c /a a ) 1/2 « 5^ » 60 A 
(£q = 11.7 A). These are the least anisotropic crystals, 
in which a crossover from incoherent to coherent inter- 
layer transport apparently takes place j|. Nevertheless, 
their anisotropy is still strongly temperature dependent 
in agreement with Eqs. (3) and (4), and changes by a 
factor of 2 - 3 between room temperature and T c . For 
other types of optimally doped layered cuprates, which 
are much more anisotropic, the condition of interlayer 
incoherence, leading to Eq. (3), is satisfied even better. 

As a second example of the applicability of Eq. (3) to 
cuprates, we consider the other doping extreme - insulat- 
ing crystals like PrBaiCu?,0-j- p0| . The incoherence of 
the c-axis transport in this case means that the localized 
states are two-dimensional (2D). In the hopping regime, 
TL is given by the average hopping distance. In the VRH 
regime uncomplicated by Coulomb interactions, the av- 
erage in-plane hopping distance R can be obtained by 
maximizing the hopping probability 
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where A is the localization length, J\f = const is the 2D 
density of states, and A a numerical coefficient. The 
only modification of the traditional treatment of VRH 
is the denominator R 2 — Rq , instead of R 2 , which takes 
into account that two localized states with close ener- 
gies cannot strongly overlap. If they do overlap, the 
phonon interaction that causes hopping will also hy- 
bridize them and push apart the energies of the new 
states. Strictly speaking, Eq. (6) is valid only for 
R Rq ~ 2A. In this limit, the average hopping dis- 
tance, which in the hopping regime determines the in- 
plane TL, is R+2R 2 /3R, where R = {XA/AfT) 1 / 3 . Then, 
according to Eq. (3), a a /a c « [4i?g/3+ (AA/A/T) 2 / 3 ]/^. 
This T-dependence of the anisotropy (a+6T~ 2 / 3 ) was ob- 
served in insulating PrBa2Cu^Oj-s and superconduct- 
ing Yi-xPrxBa^Cu^Oj-s |hJ. Note that with increas- 
ing localization length A (on approaching the MIT), the 
value of the anisotropy increases while its temperature 
dependence weakens [b/a oc A -4 / 3 ]. 

In the previous two examples we considered systems 
far away from the metal-insulator transition. The local- 
ization length on the insulating side and the correlation 
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FIG. 1. Sketch of the trajectories k(<t b ) vs. lner a . Bold 
segments correspond to "experimentally accessible" values of 
conductivities as explained in the text. The thin lines are 
hypothetical extensions of the trajectories drawn under the 
assumption that the segments corresponding to different dop- 
ing levels are parts of a continuous curve k(ct/<x), shifted with 
respect to each other due to different normalization constants 
a. The curves as shown do not represent any particular set of 
measurements, but rather a very simplified view of the overall 
behavior transpiring from available data. 

length on the metallic side of the MIT increase on ap- 
proaching the transition, and therefore, at low tempera- 
tures, the phase coherence length near MIT is dominated 
by the divergent localization/correlation length, respec- 
tively. As a result, according to Eq. (3), at a given 
T, the anisotropy p c j p a b changes nonmonotonically with 
doping, reaching a maximum in the sample closest to 
the MIT. Such a nonmonotonic variation of p c / p a b with 
doping x is observed in Y\- x Pr x BaiCuzC>T-& crystals 
p3| . At the same time, within a given range of temper- 
ature, the T-dependence of TL (and, hence, anisotropy) 
near MIT becomes substantially weaker than in either 
metallic [p c /pa ~ ^/T from Eqs. (3) and (4)] or strongly 
insulating (p c / p a ~ T~ 2 / 3 ) phase. This explains why 
the crystals of Rcfs. (2]||], which are strongly underdoped 
and much closer to the MIT than those of Ref. Jl2| , ex- 
hibit a much weaker T-dependence of resistivities and 
anisotropy at low temperatures. 

To extend the quantitative description of conductiv- 
ity to all levels of doping, we use the following idea: the 
isotropic conductance g [ Eq. (1)] of the hyperblock with 
the sides equal to the respective TL {£,£,£ } can be de- 
scribed as a function of £ only, so that the temperature 
and magnetic field dependence of g appears exclusively 
through that of £. The justification of this hypothesis can 
be partly found in the ideas of the scaling theory HQ. 
Indeed, it is known p4| that the conductivity of a conven- 



tional 3D Fermi liquid, as well as the quantum corrections 
can be expressed as functions of TL with no explicit 
temperature dependence. For a medium "composed" of 
hyperblocks {£, £, £q}, the in-plane conductivity a a differs 
from conductance g by a constant; i.e., a a (£) = g(£)/£o. 
Therefore, below we discuss the a a (£) dependence instead 
of g(£). This assumption that the T-dependence of the 
conductivities enters only through that of £, immediately 
explains why metallic o~ a and nonmetallic a c can coexist. 
For example, if a a oc f with < v < 2, then, accord- 
ing to Eq. (3), a a is metallic (da a /d£ > 0) while a c is 
nonmetallic (da c /d£ < 0). 

A quantitative description of the o~ a {£) dependence can 
be given through the logarithmic derivative 

d In a a 

K = ~dha 

which we treat as a function of cr a . This definition for- 
mally resembles the "trajectories" of the scaling the- 
ory Q. The difference is that Ohm's law does not im- 
pose any restrictions on the o~ a {€) dependence and the 
limiting values of n. Instead, we rely on experimen- 
tal data to infer the shape of the trajectory k(ct ). We 
can do this because, according to Eq.(3), the measured 
anisotropy 77 = (pc/pa) 1 ^ 2 gives directly the tempera- 
ture dependence of TL; i.e., £ = rj£ . Therefore, the 
experimental dependence a a vs. 77 determines a a (£) and 
K (f a ) = G?lncr a /dm?7 dependences. The thick segments 
in Fig. 1 are schematic representations of K(a a ) for sev- 
eral levels of doping, where the range of a a corresponds 
to "accessible" temperatures T m j n < T < T max (T max is 
typically 300 — 350 K and T m i n « T c in superconducting 
crystals) . The thin lines indicate the hypothetical exten- 
sions of the trajectories outside of this range. According 
to Eq. (3), 

d\na c „ 

— k — 2. 

d\n£ 

Therefore, both o~ a and a c are metallic for k > 2 and 
nonmetallic for k < 0, while metallic a a and nonmetallic 
a c coexist for < k < 2. 

Segment 1 in Fig. 1 represents optimally or overdoped 
regimes where both conductivities a a and a c are metallic 
at all temperatures T > T m i n because the whole segment 
is located above the threshold k = 2. The previous ex- 
ample of optimally doped cuprates [p a — a a T, £ given by 
Eq. (4), and p c — f3 c + a c T] leads to a a — q(£ 2 — £ 2 ) and 
k(ct) = 2 + a/cr, where q — ks/oiahD and a = 2q£ 2 . 

Segment 2 corresponds to a slightly underdoped sys- 
tem. At high temperatures (small o~ a ) k > 2 and, there- 
fore, <j c is metallic. With increasing £ (decreasing T) a c 
increases, reaches a maximum when n{a a ) = 2, and then 
decreases. Therefore, p c changes from metallic at high 
T to nonmetallic at low T, similar to the T-dependence 
given by Eq. (5). Within the same range of temperature, 
p a is metallic. 
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Segment 3 represents a moderately underdoped sys- 
tem. It lies entirely within the range < n < 2 and 
corresponds to metallic p a and nonmetallic p c for all 
-^miti T <i T max . 

Segment 4 corresponds to strongly underdoped crystals 
p0| , ^3[ with a a changing from metallic at high T (k > 0) 
to nonmetallic at lower T (k < 0). The singularity n = 
is integrable: 



k(<t)«±C fln^ 



1/2 



so that a a (£) as determined by the equation 



rflncr , £ 
— ; — r = In — 



(7) 



(8) 



reaches the maximum value a rnax at a finite I = 
and decreases with further increasing £ (decreas- 
ing T). Equations (7) and (8) give a a (£) « 
CTma. exp{-C 2 ln 2 (^!)/4}. 

Finally, segment 5 corresponds to an insulating crystal. 
In this case, £ f» R and cr a ~ (To exp{— 3£/A}. Therefore, 
n(a) Ri ln(cr/cr ). 

As shown, all five curves in Fig. 1 represent the same 
K,(a/a) dependence, shifted with respect to each other 
due to different values of a (determined by the density 
of carriers n and other parameters which cannot be ab- 
sorbed into TL). The reduction of n by doping reduces 
a, shifting the respective segment of k(ct/ct) to lower ab- 
solute values of a and, at the same time, ahead along the 
trajectory in terms of the reduced variable a/a. Integra- 
tion similar to Eq. (8) gives a(£) of the form 



O'a _ , 

a ~ 1 \£ 



(9) 



where f(y) is a universal function [/(l) = 1], so that 
the conductivity of the crystals at different levels of dop- 
ing is described by the same f(£) dependence, provided 
the appropriate choice of doping dependent parameters 
a and £. This hypothesis can be verified experimentally 
by checking whether the conductivity at different dop- 
ing levels plotted against anisotropy rj form a continuous 
curve when a a and rj are properly normalized. 

This idea of universality of the a a {£) dependence allows 
the prediction of the behavior of the normal state conduc- 
tivities at temperatures below T c . For example, let seg- 
ment 3 in Fig. 1 represent a superconducting crystal with 
metallic a a and nonmetallic a c for all T c < T < T max . If 
we suppress the onset of superconductivity with a mag- 
netic field in order to reveal the normal state at T < T c , 
we will observe that the trajectory crosses the thresh- 
old k = and a a also becomes nonmetallic. This type 
of development has been observed |^],[| in underdoped 
cuprates. 

It is possible that the trajectories shown in Fig. 1 
oversimplify the situation, and there may be more than 



one type of f(£) dependence. The data suggest that in 
Yx-xPi'xBavCusO'TS, a metal-insulator transition takes 
place at intermediate doping levels somewhere between 
optimally doped and strongly insulating phases As 
a result, there may be at least two types of trajectories 
(or a branching point), corresponding to the phases on 
either side of the MIT. 

The assumption that the response of the conductiv- 
ities to external perturbations other than temperature 
also comes from that of TL allows us to predict these 
responses because 



dlna a dint <91ner c . „. dint 
k — — — ; - - = (k-2)- 
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dx 



dx 
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Here x can be magnetic field or concentration of impu- 
rities. For k > 2 or n < both conductivities change 
similarly. However, for a wide range of doping levels and 
temperatures [Fig. 1], for which < n < 2, the two 
conductivities respond oppositely. 

For example, in Bi2Sr2Ca(Cui- x Zn x )20s+ y single 
crystals a c increases and a a decreases with increasing 
Zn concentration ]l5[ ]. The impurities reduce the in- 
plane TL (d£/dx < 0) by decreasing the elastic mean 
free path and the diffusion coefficient. The crystals used 
in |]l5| exhibit metallic p a (T) and nonmetallic p c (T), so 
that < k < 2 and, as a result, da a /dx < while 
da c /dx > 0. Sometimes, the increase of the c-axis con- 
ductivity in response to a perturbation has been inter- 
preted in literature as a crossover to coherent transport 
in the c-direction. We see that this is not necessarily the 
case. The anisotropy may decrease and a c increase due 
to the reduction of the in-plane TL, even when the c-axis 
TL remains fixed. 

The effect of the magnetic field H on conductivities is 
determined by the destructive interference Q, resulting 
in the reduction of TL (d£/dH < 0). Therefore, accord- 
ing to Eq. (3), the magnetoanisotropy [S(p c /p a )] should 
be negative, while the sign and relative magnitude of the 
magnetoconductivities Aa a and Aa c depend on the value 
of k according to Eq. (10). One can see that the magne- 
toresistivities correlate with the temperature coefficient 
of the respective component of the resistivity: 



dp a ,c 
OH 



= Q 



dpa :C 

dT 



where Q = (d£/dH)/(d£/dT). Since Q > 0, the sign 
of magnetoresistance Ap a c is the same as the sign of 
dp a ,c/dT, as indeed was reported in Refs. ||,§,|l6|. The 
magnitude of the magnetoeffects depends on the number 
of flux quanta (hi) permeating the hyperblock. For H\\c, 
h c = H£ 2 /(f>o = ( Pc /pa)H/H (H = (f> /£l and O is the 
quantum of flux). If the field is parallel to the planes, 
h a = H££ /<j> = {p e /Pa) 1/2 R/H . Thus, h c » h a which 
may result in a strong angular dependence of the mag- 
netoeffects. 
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